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ABSTRACT 

We estimate lepton capture and emission rates, as well as neutrino energy loss 
rates, for nuclei in the mass range A=65-80. These rates are calculated on a tem¬ 
perature/density grid appropriate for a wide range of astrophysical applications 
including simulations of late time stellar evolution and x-ray bursts. The basic 
inputs in our single particle and empirically inspired model are i)experimentally 
measured level information, weak transition matrix elements, and lifetimes, ii) 
estimates of matrix elements for allowed experimentally-unmeasured transitions 
based on the systematics of experimentally observed allowed transitions, and 
iii) estimates of the centroids of the GT resonances motivated by shell model 
calculations in the fp shell as well as by (n,p) and (p,n) experiments. Fermi 
resonances (isobaric analog states) are also included, and it is shown that Fermi 
transitions dominate the rates for most interesting proton rich nuclei for which 
an experimentally-determined ground state lifetime is unavailable. For the pur¬ 
poses of comparing our results with more detailed shell model based calculations 
we also calculate weak rates for nuclei in the mass range A=60-65 for which 
Langanke and Martinez-Pinedo have provided rates. The typical deviation in 
the electron capture and beta decay rates for these ~ 30 nuclei is less than a 
factor of two or three for a wide range of temperature and density appropriate 
for pre-supernova stellar evolution. We also discuss some subtleties associated 
with the partition functions used in calculations of stellar weak rates and show 
that the proper treatment of the partition functions is essential for estimating 
high temperature beta decay rates. In particular, we show that partition func¬ 
tions based on un-converged Lanczos calculations can result in estimates of high 
temperature beta decay rates that are systematically low. 


1. Introduction 


In this paper we provide estimates for weak interaction rates involving intermediate 
mass nnclei. Anfderheide et al. (1990) and Anfderheide et al. (1994) have argued that at 
late times the electron fraction in the Fe core of pre-supernova stars can be so low that weak 
processes involving the A> 65 nuclei we study are important. Depending on the entropy per 
baryon, which determines the free proton fraction, electron capture on heavy nuclei may also 
play an important role during collapse (Bethe et al. 1979; Fuller 1982). In addition, the weak 
rates we provide for proton rich nuclei may be used in studies of nucleosynthesis and energy 
generation in x-ray bursts and other rp-process sites (Wallace and Woosley 1981). For the 
rp-process, weak rates are needed for proton rich nuclei at least up to mass 110. Electron 
capture and positron decay rates for proton rich nuclei in the mass range A=81-110, as well 
as a discussion of some peculiarities of the weak rates in the rp-process environment, will be 
presented in another paper (Pruet and Fuller, in preparation). 

The formidable task of calculating the electron and positron capture and emission rates 
in the conditions characteristic of these astrophysical environments has received more than 
four decades of attention. The hrst self consistent calculations to include the effects of 
the Fermi and Gamow-Teller (GT) resonances as well as the thermal population of these 
resonances for a broad range of nuclei and thermodynamic conditons were done by Fuller, 
Fowler and Newman (1980, 1982,b) (hereafter FFN), Fuller (1982), and Fuller, Fowler and 
Newman (1985). FFN presented a physically intuitive and computationally tractable method 
for determining the strength and excitation energy of the Fermi and GT resonances. The 
groundwork for the treatment of these resonances in a thermal environment was laid down 
by Bethe et al. (1979) and FFNII. The importance of hrst forbidden transitions, blocking, 
and thermal unblocking for the very neutron rich nuclei present during collapse was pointed 
out by Fuller (1982). Gooperstein and Wambach (1984) calculated electron capture rates 
for these nuclei by estimating the parity forbidden matrix elements as well as the effects of 
thermal unblocking of the allowed GT"*" strength. 

The last two decades have seen a great increase in our understanding of weak inter¬ 
action systematics in intermediate mass nuclei. There are now semi-direct measurements 
of the GT strength distribution for roughly 40 nuclei from forward angle (n,p) and (p,n) 
scattering experiments. There are also shell model calculations for the strength distribu¬ 
tion for some 400 nuclei in the fp shell (Gaurier et al. 1999; Langanke and Martinez-Pinedo 
2000), as well as a number of RPA and QRPA calculations of the GT resonance in heavier 
nuclei. One consequence of these studies is that the FFN prescription misses some impor¬ 
tant systematics in the assigning the centroid of the GT resonance. The influence of this 
misassignment on the rates was discussed by Aufderhide et al. (1996), Gaurier et al. (1999), 
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and Langanke and Martinez-Pinedo (2000) (hereafter LMP). Langanke and Martinez-Pinedo 
(2000) also provided npdated calcnlations of the rates for A< 65 based on large scale shell 
model calculations. 

We attempt to remedy the misassignment of the GT centroids. Other than this and 
a different treatment of high temperature partition functions, our strategy for calculating 
weak rates for A>65 is essentially the FFN approach. In this approach the rates are broken 
into two pieces: a low part consisting of discrete transitions between individual levels, and 
a high part involving the Fermi and GT resonances. This approach is valid provided that 
discrete transitions between high lying levels never dominant the rates. In addition, available 
experimental information for log ft values and level energies, parities, and spins is used. The 
inclusion of this much data for more than a hundred nuclei is a difficult task, and would not 
be possible without the web-based nuclear structure databases (in particular those provided 
by nudat and the table of isotopes). 

The simple FFN approach for estimating weak rates may be the most natural framework 
for incorporating experimentally-determined nuclear properties. Additionally, a virtue of 
using this semi-empirical schematic approach is that we can easily see where the key nuclear 
uncertainties are. Our work can then serve as a catalyst for more detailed follow up nuclear 
structure studies for important rates. 

In the next section we present a discussion of the formalism for calculating weak rates. 
In section III we discuss the assignment of experimentally unknown matrix elements for 
allowed, discrete state transitions. This discussion is based on the systematics of experimen¬ 
tally observed weak transitions. In section IV we present a simple approach for estimating 
the position of the GT resonances and make some comparison with data and shell model 
calculations. The calculation of rates in a high temperature environment is treated in section 
V. We also discuss in section V some subtleties associated with the proper partition function 
to be used in these calculations. Section VI gives some comparisons of our estimates for the 
rates with those provided by Langanke and Martinez-Pinedo. Finally, we conclude with a 
discussion of the results. 


2. Weak rates formalism 

The total decay rate for a nucleus in thermal equilibrium at temperature T is given by 
a sum over initial parent states i and hnal parent states j; 

A = Pi Xij , 

i 3 


( 1 ) 
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where the population factor for a parent state i with excitation energy Ei and angular 
momentum Jj is 

, 2 ) 

z 

with 


Z — + 1)( 


-Ei/kT 


the nuclear partition function. Here Xij is the specific weak transition rate between initial 
parent state i and daughter state j and is formally given by 

~ Wur ^ ^ 


where 


ifth m, {ftw 


and where the relation between ft-values and the appropriate Gamow-Teller or Fermi matrix 
element is 


Here the total Gamow-Teller matrix element between initial parent state {ipf) and hnal 
daughter state I'lpf) is 

\Mgt\1 = I 5^c^n(r±)n|^f)r (8) 

n 

where the sum is over all nucleons. The GT strength satishes the sum rule 

V - Sp+ = - | (^/1 ^ a„(r+)„|^f) 1^ = 3{N-Z). (9) 

f n n 

Similarily, the Fermi matrix element is 

|a/fP = 15^(T±)„|i/>f )l" = t(t +1) - r,(r. -i) = \n-z\. (lo) 

n 

This equation is derived by noting that [T+,T_] = 2Tz = {N — Z). 

The phase space factors for {3^ decay are 


fij= / w‘^{qn-wfG{±Z,w){l-f^){l-f^,)dw 
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and for electron/positron capture are 


fij 


w‘^{qn + wfG{±Z, w)/=f(1 - fu)dw 


• Wl 


( 12 ) 


where the upper(lower) signs are for electrons(positrons), = (Mp — + -E* — Ej)/meC^, 

and where Mp is the nuclear mass of the parent and is the nuclear mass of the daughter. 
The threshold is wi = 1 for g„ > — 1 and wi = |g„| for < —1. The lepton occupation 
factors are 


h 


U-U^ 



(13) 


where the (-) sign is for electrons, the (+) sign is for positrons, U is the electron kinetic 
energy and Up the kinetic chemical potential (total chemical potential is here dehned to 
include the electron rest mass). The factor G{±Z,w) is the coulomb wave correction factor 
dehned in terms of the Fermi factor F 


G{±Z,w) = -F(±Z,w) 
w 

as discussed in FFNI. 


(14) 


3. Assignment of unmeasured, allowed transition matrix elements 

As in FFN, we break the rates Xij into two distinct components: non-resonant discrete 
transitions between low lying levels, and transitions involving resonance states carrying total 
weak strength (|MgtP + I^fP) of order one or greater. Discrete state transitions between 
low lying states in the parent and daughter are important when the capturing lepton energies 
are too low to reach resonance states in the daughter, or when the temperature is too low to 
thermally populate parent levels with fast transitions to daughter states. To estimate matrix 
elements between low lying parent and daughter levels, FFN adopted the simple prescription 
that all transitions not forbidden by the selection rules have some average matrix element 
characteristic of a group of nuclei. The validity of this approach can be addressed by looking 
at the systematics of matrix elements for experimentally observed j3^/ec and (5~ decays. 

Figures 1, 2, and 3 show the characteristic log ft value for all nuclei in the mass range 
A=65-80 for which there is experimentally-determined weak decay information. By ” charac¬ 
teristic log ft" we mean here the log ft value obtained by assuming that all of the measured 
discrete strength is spread uniformly over all states in the daughter for which the selection 
rules do not forbid a transition and for which the Q-value is positive. This estimated number 
of allowed transitions (uaiiow) for each nucleus is shown on the x axis. The determination 
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of rzaiiow is difficult because of uncertainties in the angular momentum and spin assignments 
of levels. Where the angular momentum of a given level is uncertain we have adopted the 
middle value if more than two possibilities are listed, and the largest value if only two pos¬ 
sibilities are given. Where the parity is listed, but uncertain, we have adopted the tentative 
value. Transitions involving a level for which no spin or parity information is given are 
labelled as forbidden. 

For illustration, consider the (3 decay of ®®Co. Experimentally it is observed that ®®Co 
decays in this channel to two states in ®®Ni, one with a log/t of 4.21 and the other with 
log/t =4.75. By examining the experimentally studied levels in ®®Ni it is seen that there 
are 7 levels that have spins and parities consistent with allowed decay from ®®Co. The 
characteristic log ft for ®®Co is then dehned as — logio((l/7)(10“^-^^-|-10“^-’^)). Our dehnition 
of the characteristic log ft is rough in that it neglects the experimental difficulties associated 
with measuring weak transitions. In particular, uncertainties arising from the possibility 
of the feeding of daughter states from higher lying states has not been accounted for, nor 
has the difficulty of observing near-threshold transitions. Nonethless, a case may be made 
from these hgures that the assumption of an average matrix element is a reasonable one, 
particularly when several transitions are involved. For example, for nuclei with 65 < A < 70, 
only one nucleus with Uaiiow > 2 has a characteristic log/t differing from 5.4 by more than 
1, and for nuclei in the range 70 < A < 80 only 5 nuclei with Uaiiow > 2 have a characteristic 
log/t differing from 5.7 by more than 1. 

In this work we take a characteristic log/t of 5.4 for nuclei with A<70, and a character¬ 
istic log/t of 5.7 for nuclei with A>70. An alternative approach for estimating these matrix 
elements would be to assign log/t’s from a statistical distribution. Our procedure should 
give a reasonable estimate of the rates when several transitions contribute nearly equally to 
the rate. However, our estimated rate is obviously subject to uncertainty when only one or 
two experimentally unknown transitions dominate. An important class of such nuclei are the 
proton-rich even-even nuclei with nearly closed neutron and proton sub-shells. These can 
have anomolously large 0^ —W GT transitions. In x-ray burst environments, with tem¬ 
peratures T ~ 0.25 MeV, the first 2^ excited state in these nuclei is thermally populated. 
The thermal j3^ decay rate of the nucleus depends sensitively on whether this 2^^ state also 
decays with anomolously large matrix elements. Schatz et al. (1998) studied this question 
for a number of proton-rich nuclei and found that the 2+ state decays more rapidly than 
the ground state when the ground state decays quickly. We incorporate the results of the 
Schatz et al. (1998) shell model study for the four nuclei (72,74j,^j,^76,78 -with fast 0+ —I’*' 
transitions. For these nuclei we assume that the first 2+ state has a strength distribution 
identical in shape, but with matrix elements twice as large, as the strength distribution from 
the O’*" ground state. This gives decay rates within ~ 30% of those calculated by Scahtz et 
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4. The Fermi and Gamow-Teller Resonances 

The Fermi resonance \Fi) corresponding to a given state {ipi) is generated by application 
of the isospin raising or lowering operator, \Fi) = T±\ipi) (Eq. 10). The selection rules for 
Fermi transitions are Att = AT = AJ = 0. Because AT = 0 and because the ground 
state of a nucleus generally has the lowest possible isospin, there typically is only non-zero 
Fermi strength for transitions from a nucleus with greater isospin T^ to a nucleus with 
lesser or equal isospin T^. Since the nuclear part of the Hamiltonian (Rnuc) is isospin 
independent and the electromagnetic part is small in comparison, the resonance generated 
by T± is narrowly concentrated about the IAS. The excitation energy can be estimated from 
the difference in Coulomb binding energy of the parent and daughter nucleus. FFNI gives a 
useful approximation for the excitation energy of the IAS in the daughter nucleus: 

Tias = Ap-Ad T 0.7824 ± 1.728 min(Zp, Zd)/R (15) 

where R « + 0.78 is the nuclear radius in fm, p and d refer to the parent and 

daughter, respectively, A is the atomic mass excess. The upper signs in the above equation 
correspond to (Z, N) —> (Z -|- 1, N — 1) transitions for neutron rich parents, while the lower 
signs correspond to (Z, N) —(Z — 1, N -|- 1) transitions for proton rich parents. Eq. 15 agrees 
well with measured and shell model predictions for IAS energies. 

The Gamow-Teller operator is 

GT± = ^o-r^(n), (16) 

n 

where the sum is over all nucleons n. The collective GT resonance state |CGTi) corresponding 
to a given parent state \'ipi) is given by application of the GT operator, |GGTi) = GT|'^j). The 
selection rules for GT transitions are Att = 0, AJ = 0, ±1 no 0 —0, and AT = 0, ±1. The 
GT strength distribution is harder to characterize because iJnuc is strongly spin dependent. 
Since [Jfnuc, GT±\ ^ 0, the GT strength can be fragmented over many daughter states. 
However, in practice the stellar weak rates are usually determined by the total strength and 
the centroid of the strength in excitation energy (provided that low lying discrete transitions 
are well accounted for). 

The strength in the GT resonance was also estimated by FFN in a zeroth order shell 
model picture. In this picture the lowest shell orbitals are hlled with nucleons, and the total 
strength is taken to be the sum of the contributions from each pair of single particle orbitals: 
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Here i and f denote initial and final orbitals respectively, Up and Uh denote the number of 
particles and holes in these orbitals, and is the single particle matrix element connecting 
the initial and final states. These single particle matrix elements can be found from angular 
momentum considerations and are shown in Table 1 of FFNII. 

We follow FFN in using the single particle result (Eq. 17) to estimate the total strength. 
Experimentally it is well established that the axial vector current is renormalized by a factor 
of ~ (1/1.24) in nuclei. This results in a strength a factor of (1/1.24)^ 1/2 smaller 

than shell model calculations give. In addition, shell model calculations show that residual 
interaction-induced particle-hole correlations further reduce the total strength by a factor of 
one to a few. Typically, these correlations are more important for GT’*' transitions, so that 
the additional quenching is larger for these transitions. We adopt a quenching factor of 4 for 
GT"*" transitions and 3 for GT“ transitions. These values for the quenching factors generally 
give strengths within a factor of two of more detailed strength determinations. When the 
quenched value for the strength is less than one, we assign a configuration mixing strength 
of one. As discussed below, this is roughly consistent with the results of (n,p) experiments 
for GT"*" blocked nuclei. 

FFN estimated the centroid of the GT resonance by considering a zeroth order shell 
model description for the spin flip part of the GT resonance. This configuration was com¬ 
pared with the zeroth order shell model description of the daughter ground state and assigned 
an excitation energy 

F/qt = AF/s,p. -|- AEpair + AEph. (18) 

Here AEg.p, is the difference in single particle energies between the two states (daughter 
ground state and spin flip GT resonance state), AEpair accounts for the difference in pairing 
energy between the two states, and AEph (taken to be 2 MeV) accounts for the effects of 
configuration mixing and particle hole repulsion. For transitions the energy of the 

IAS in the daughter is added to eq. 18. 

For the assignment of the GT centroids we adopt a procedure close to that outlined 
by FFN. However, as noted in the introduction, the FFN approach misses some important 
systematics in the centroids of the strengths as revealed by more recent experimental data and 
shell model calculations. One potential remedy for this is to do RPA or QRPA calculations 
for the strength distributions for nuclei too heavy to be studied via detailed shell model 
calculations. Such calculations have been done for a number of nuclei by several different 



9 


groups. A simpler approach is to approximately account for the effects of the competition 
between the k^ri^ri ■ Ti)(a 2 ■ T 2 ) and kr{Ti ■ T 2 ) terms in the nuclear Hamiltonian. The effect 
of this competition on the centroid of the GT resonance is most easily seen in the Tamm 
Dancoff approximation from the argument given by Bertsch and Esbensen (1987). In this 
approximation the kr and k^r forces give rise to a GT“ excitation energy scaling as 

A'gT- ~ A-IAS + ^Es.o + ‘^{karScT- — {N — Z)kr)- (19) 


Here AEg.o. is the spin orbit splitting characteristic of the single particle transitions for 
the GT“ resonance. As an example of this equation, consider the case of ^"^Fe and ^®Fe. The 
zeroth order (3~ strengths for these nuclei are |MgtP = 16.3 ( ^^Fe) and 22.3 (^®Fe). For 
k^^ ~ (20/A)MeV ~ (2/3)A;^, Fqn. 19 implies {Egt- - FiAs)56Fe - (^gt- - FiAs)54Fe ~ 
— 1 MeV. Shell model calculations and (n,p) experiments for these nuclei give this difference 
as approximately -2 MeV. As Fe becomes more and more neutron rich Fqt- approaches Fjas 
and eventually falls below it (although perhaps beyond the neutron drip line for a nucleus 
as light as Fe). 

A somewhat more sophisticated approach to incorporating the effects of the competition 
between the kr and k^r forces is the random phase approximation with a separable force. In 
this approximation the GT’*' and GT“ resonances become eigenstates of the Hamiltonian. 
The energies of the resonances are approximately given by the roots of the algebraic equation 


fl ^ J2_ ^ h ^ 

e* — e + A“q Q — e —Ci + e + A+q 2 /cg-r *S'tot 


( 20 ) 


(see Gaarde et al. (1981) for an application of this equation to experimentally observed 
P~ strengths, or Rowe (1970) for a more pedagogical discussion). In Eq. 20, Atotai = 
>S'gt- + *S'gt+ = 3\N — Z\ + 2S'gt+ is the sum of the strengths in the plus and minus 
directions, /i is the fraction of this strength in the GT“ spin-flip mode, /2 is the fraction of 
this strength in the GT“ non-spin flip mode, and /a is the fraction of strength in the GT"*" 
direction. The spin orbit splittings A+q and A“q are the splittings appropriate for the spin 
flip transitions in the plus and minus directions, e* = (e^ — ^u)i is the difference in the proton 
and neutron single particle energies for the levels involved in the transition. The quantity 
Ci is related to the energy of the IAS by Fias — = 2/ct-|V — Z\. The largest root of Eq. 

20 corrresponds to the energy of the spin flip mode. Eq. 20 reduces to Eq. 19 in the limit 
/i = l. 


Our approach for calculating the centroid of the GT“ resonance is based on Eq. 20. 
The strengths in this equation are estimated in the zeroth order shell model picture de¬ 
scribed above. For consistency with FFN we take the spin orbit splittings from Seeger and 
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Howard (see the table in Hillman and Grover (1969)). When more than one spin flip tran¬ 
sition contributes to the strength we take the strength-weighted average. The parameter k^ 
can be estimated using measured IAS energies and estimates for the particle-hole energies. 
Alternatively, an estimate for the IAS energy (Eq. 15) gives a relation between kr and e*. 
The parameter k^r can be chosen to give good agreement with shell model and experimental 
results for the GT“ resonance in the fp shell. In this work we adopt kr = 28.5/A MeV and 
krrr = 23/A MeV. These values are close to those given in Bertsch and Esbensen (1987) and 
Gaarde et al. (1981). It has previously been noted by a number of authors (e.g. Gaarde et 
ah (1981)) that a simple prescription can do a fair job of predicting the centroids of the GT“ 
resonances. In table 2 we reaffirm this by comparing the predictions based on Eq. 20 with 
measured and shell model results. 

For the GT^ resonance, estimates based on a separable force are not well justihed. 
Higher order particle hole correlations and correlations induced by other terms in the Hamil¬ 
tonian play a more important role. Nonethless, the RPA result accounts for the spin orbit 
splittings and the systematics of the influence of the kr and k^r forces in an approximate 
way. We estimate the centroid of the GT"*" resonance from the equation analogous to Eq. 
20 (which is found by reversing the role of the plus and minus transitions, and by setting 
Ci = — Ctt) or alternatively by setting e ——e in Eq. 20). We also add an additional term 

to account for the effects of correlations missed by the separable force estimate: 

EgT+ = -SgT+,RPA + <^ph- (21) 

A value of 2 MeV for the empirical correction (5ph gives good agreement with the results 
of experimental and shell model studies of nuclei in the lower half of the fp shell. This is 
demonstrated in Table 3. 

Our simple estimate for Egt+ probably breaks down as the fp shell approaches being 
hlled. Fortunately, for most nuclei in this case {A > 65) the GT"*" transition is nearly blocked 
(see below). Likewise, for those proton rich nuclei with substantial amounts of unblocked 
strength the electron capture Q-value is large, so again the rate is not terribly sensitive to 
assumptions about the centroid. For the few proton-rich nuclei with modest GT"*" strength, 
we can compare our results to the more detailed QRPA calculations of Sarriguren et al. 
(2001). These authors performed QRPA studies for even-even Ge, Se, Kr, and Sr isotopes. 
For the most part their calculated electron capture half lives agree well with experimental 
half lives. 

For ®®Ge, ®®Ge, and ^'^Se Sarriguren et al. (2001) hud that the strength distribution is 
not very sensitive to which of two nearly degenerate shapes the nucleus assumes, so that 
a direct comparison with our estimate is sensible. For these nuclei our estimate for the 
centroid of the resonance is lower than theirs by 1.5 Mev(®®Ge), 1.8 MeV (®®Ge), and 2.5 
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MeV (^°Se). For some of these nuclei it is not clear which estimate is correct. Sarriguren et 
ah (2001) show that their method gives a centroid about 2 MeV higher than the experimental 
centroid for ®^Fe and ®®Fe, and that there calculation for ^°Ge misses a modest amount of 
experimentally-determined low lying strength. 

However, in at least one case C^^Se), our assignment of the strength is at least 0.5-1 
MeV too low, because it results in too much strength within the electron capture window 
and results in a half life about 5 times shorter than the experimentally determined half life. 
(Note though, that where a half life is available for a nucleus our calculations always agree 
with that half life. We do not include resonance strength below the Q-value window for a 
transition when experimental information is available.) We will argue in the last section that 
the uncertainty in the placement of the centroid of the GT’*' strength is not very important 
for nuclei in the mass range we are considering. 

As the GT+ strength decreases (i.e., as an isotope becomes more neutron rich), Eq. 
21 eventually gives a negative excitation energy for the resonance in the daughter. This 
typically happens when the single particle estimate for the strength (Eq. 17) is less than 
about 5. In this case, the strength is dominated almost entirely by conhguration mixing. 
There are a few (n,p) studies of such nearly blocked nuclei. Vetterli et ah (1992) studied 
^°Ge(n, p)^'^As and ^^Ge(n, p)'^^As. The experimentally-determined GT strength for ^°Ge is 
B(GT)=0.84 ± 0.13 or B(GT)=0.72 ± 0.14. The two different values correspond to different 
ways of estimating the AL = 0 component of the (n,p) cross section. The higher estimate 
comes from a multipole decomposition (m-d) of the cross section, while the smaller estimate 
is derived by approximating the cross section measured at 5.8° as the AL = 1 component 
of the cross section. For the m-d the strength distribution for ^‘^Ge is approximately flat up 
to a few tens of MeV in excitation energy in ™As. For the 5.8° subtraction method, the 
strength falls after about 6 MeV in excitation energy in ^'^As. For ^^Ge B(GT)=0.23 ± 0.05 
(5.8°) or B(GT)=0.86 ± 0.14 (m-d) and the strength distribution is roughly flat or falls off 
after about 6 MeV depending on the background substruction method. Helmer et ah (1997) 
have studied ^®Se(n, p)^^As. They hnd B(GT)=1.45 with a roughly flat strength distribution 
for the m-d method, and B(GT)=0.35 with a strength extending to about 6 MeV in ^®As if 
the 6° data is used to estimate the AL = 1 component of the cross section. In this work we 
adopt the procedure that when the RPA estimate Eqx+,rpa appearing in Eqn. 21 is negative 
with respect to the daughter ground state, the strength distribution is represented by a 1 
MeV gaussian centered at 1.8 MeV. This is approximately consistent with the experimental 
data analyzed using the 5.8° (6°) background substruction method. Our prescription misses 
the high lying strength estimated from the experimental data analyzed using the multipole 
decomposition method. We will discuss in the last section the uncertainties in the rates 
resulting from the unknown strength distribution for these blocked nuclei. 
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5. Thermal considerations 

The picture for Fermi and Gamow-Teller strengths outlined above becomes more com¬ 
plicated at high temperatures on account of the thermal population of parent excited states. 
In evaluating the contribution to the rates of transitions between low lying levels, we use the 
same set of experimentally-determined levels discussed above. 

As the temperature rises, the evaluation of Eq. 1 requires knowing the strength dis¬ 
tributions of an impossibly large number of states. For example, at temperature of 1 MeV 
the partition function for an odd-odd nucleus with A 60 is a few hundred, and the mean 
excitation energy is ~ T‘^A/8 ~ 7 MeV (here pp ^ a is the level density at the fermi 
surface and a ~ A/8 is the level density parameter). The approximation traditionally (see 
FFNII) used to make the problem tractable is the Brink approximation, which postulates 
that the centroid of the Gamow-Teller strength distribution corresponding to a parent state 
at excitation energy Ei is shifted up by an energy Ei with respect to the centroid of the 
strength distribution corresponding to the parent ground state. It is generally assumed that 
the total strength remains the same for all transitions. The validity of the Brink approxi¬ 
mation has been investigated in some detail by BMP. They hnd that the approximation is 
good for the first few low lying states for which they calculate strength distributions. 

With the Brink approximation, the contribution of discrete state to high lying resonance 
state transitions can be approximated. For dehniteness we assnme in the following discus¬ 
sion that the parent nncleus has isospin T>, the daughter nucleus has isospin T<, and that 
the GT“ operator acting on the parent generates states in the daughter (and conversely). 
Each parent state (|'^f)) has a corresponding Fermi resonance (IF)) = T~\iljf)) and collec¬ 
tive Gamow-Teller resonance {\CGTi) = GT~\ipf)) in the daughter. With the Brink 
approximation, the Q-valnes for these transitions are independent of the excitation energy of 
the thermally populated parent state. The overall transition rate for these processes can be 
obtained brom a calculation of the gronnd state rate alone, but with the population factor 
of the gronnd state set to nnity. This is part of the “FFN trick”. 

Likewise, it is possible to thermally populate the collective Gamow-Teller states in 
the parent {\GGTj) = GT~^\'ipf)). These thermally populated states can decay to the 
danghter with large overlap. These “back-resonances” are shown in Fig. 4. Suppose for a 
moment that the GT"*" resonance corresponding to a state in the daughter is conhned to a 
single state in the parent. In this case the Brink approximation implies that for each state 
with excitation energy Edaughterj in the danghter, there is a corresponding GT+ resonance 
state {\GGTj) = GT+|'^°)) in the parent with excitation energy EoS-Edaughterj- Here Eq is the 
excitation energy in the parent of the resonance state corresponding to the danghter gronnd 
state. Because in the Brink approximation the Q-valnes and strengths for all the transitions 
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{Eq + -Ej)parent “> (-^j)daughter a^e the Same, the rate contribution for each transition is also 
the same and can be written as A. The contribution of the back resonances to the rates can 
be directly evaluated: 


A 


backres — 


^ 6Xp( -Eparentii/^) (2'^parent,i “1“ 1) 


-^parent 


( 22 ) 


^ j ( ( -EdaughterJ -E'o)/^) (2(/daughter,] H" 1) 


■^parent 


= Xexp{-Eo/T) 


■^daughter 

7 

•^parent 


Eq. 22 was first derived in FFNII and is sometimes referred to as the “FFN trick”. 


In this work we represent the resonance strength distribution by a gaussian of nominal 
width a. For the GT’*' strength distribution we take a = 1 MeV and for the GT“ strength 
distribution we take a = 2 MeV. This approximation does not fairly represent the complex 
structure and variation seen in real strength distributions. However, we will see in the last 
section that this simple approximation captures the features relevant for calculations of weak 
rates in a stellar environment. The case where the resonance is spread over several states 
in the parent can be treated similarity to the case where the resonance is concentrated in a 
single state. If the resonance corresponding to the ground state in the daughter is spread 
over states with excitation energy Eq, Ei, ...En in the parent, with corresponding rates to 
the daughter ground state given by Aoj, then Eq. 22 becomes 


A 


AA,- Aq.) exp( F/gj/T) .^(jaughter 


backrate 


-^parent 


(23) 


This equation is exact in the limit where the Brink approximation holds (i.e., the strength 
distribution from every daughter state is identical to the strength distribution from the 
daughter ground state except for an overall shift in energy). The interesting feature of Eq. 
23 is that if all of the Q-values for the n states are comparable, then 


Abackrate ~ ^(Aoj) exp( F/g/T).^i-[aughter/-^parent Aexp( F/g/F').^i-[aughter/-^parent) (2^) 

independent of n. Here the relation n{\oj) ~ A (where A is the rate of decay calculated 
assuming that the GT resonances are conhned to single states) holds because the total 
strength in the GT resonance is independent of the number of states it is spread over. In 
a sense Eq. 24 is counterintuitive because one might expect that if the strength is spread 
out over n states in the parent, the decay rate from the parent should be a factor of 1/n 
smaller. However, in some cases the strength distributions arising from different states in 



the daughter must overlap (i.e. correspond to identical states in the parent) to the extent 
that the level densities in the parent and in the daughter are comparable. This is accounted 
for by the ratio of partition functions appearing in Eq. 23. 

In principal, the back resonance contribution to the rates could be calculated through 
a direct (e.g., shell model) calculation of the strength distribution from the daughter states. 
In practice this approach is not feasible because such a calculation would have to have fully 
converged final states corresponding to the strength distribution from many daughter states 
in order to distinguish which ones are identical. Otherwise, an overcounting of the partition 
function results, and the hnal rate estimate is a function of how many Lanczos iterations are 
done. This line of reasoning indicates that at high temperatures the current shell model-based 
calculations could systematically underestimate the contribution to the (3 decay rates from 
the decay of the thermally populated GT+ resonance states by up to an order of magnitude. 
To show this, we plot in Fig. 5 a comparison between our calculated f3~ decay rates and those 
of LMP at a temperature of Tg = 30 and a density of pYe = lOg/cm^. These thermodynamic 
conditions are artihcial, but serve the purpose of illustration. When the f3 decay rate is fast 
(and insensitive to the hner details of the strength distribution), the discrepancy between 
the two calculations is about a factor of 10. This is roughly the number of independent 
Lanczos states carrying strength per daughter state in the calculations of LMP. From our 
simple analysis it isn’t clear that either set of rates is more reliable. However, it is clear 
that differences in treatments of the partition functions can result in signihcant differences 
in estimates of the rates. 

The trend of greater disagreement with faster decay rate seen in Fig. 5 arises from the 
competition of two factors. For low decay rates, the Q-value of the decay (excitation energy 
of the GT+ resonance with respect to the daughter ground state), is typically small. In 
this case the width of the resonance enhances the decay rate compared to the decay rate 
calculated from our artihcially narrow resonance. For faster decay rates, the Q-value is large, 
the lifetime is relatively insensitive to the width of the resonance, and the Lanczos-based 
calculation rate estimate is simply suppressed by the overestimate of the parent partition 
function. 

The effect, and the importance of consistency in partition functions in general, can 
be illustrated by considering the conditions in the post-silicon burning, pre-collapse core 
of a massive star. There electron capture proceeds on iron peak nuclei, driving them to a 
neutron excess where “reverse” (3~ decay balances “forward” electron capture. For example, 
Aufderheide et al. (1994) identify ®^Gr as the endpoint nucleus where the forward and reverse 
neutronization rates balance. In balanced conditions, partition functions are crucial rate and 
abundance determinants. We can compare our rates with LMP as in Fig. 5, but now for 
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Tg = 10, pYe = 10®g/cm^, roughly approximating immediately pre-collapse conditions. This 
comparison is shown in Fig. 6. Again there is a systematic trend: the LMP rates are lower 
than ours by a factor of 4 on average with a fair scatter. This is smaller than the disagreement 
presented above for the fastest f3 decay rates, but still of potential signihcance given the 
dependence of the initial Fe core mass on the electron fraction. Part of the discrepancy 
between the rate estimates undoubtedly stems from differences in placement and width of 
GT strength and other nuclear uncertainties. 

However, the partition effect outlined above likely plays a role as well. In fact, these 
conitions (Tg = 10, pYe = 10®g/cm^) are electron degenerate with Fermi energies ~ 5 MeV, 
precluding signihcant (3 decay for all but the nuclear decay pairs with relatively large Q- 
values, just those where we argued that the partition function-based uncertainties in rate 
estimates could be large. 

In evaluating the ratio .^daughter/-^parent in Eq. 23 we use the compilation of parti¬ 
tion functions from Rauscher and Thielemann (1999). These partition functions include 
experimentally-determined low lying levels and are supplemented at higher excitation ener¬ 
gies by a level density calculated from a back shifted Fermi gas formula. At low temperatures 
the partition functions of Rauscher and Thielemann (1999) agree well with the partition func¬ 
tions we calculate for evaluating the rates between low lying levels. For temperatures above 
2 MeV we take .^daughter/-^parent = 1- This is valid because the mean excitation energy at 
these temperatures is well above the pairing gap. 

We do not claim that our partition function treatment is necessarily better than others, 
and it may well be inadequate in some conditions. In fact there is a basic inconsistency 
in our rates: we include many states in our partition function sums for which we include 
no weak interaction strength. Ideally we should include all states and all associated weak 
strength: only fully converged Lanczos and Monte Carlo calculations of weak strength and 
partition functions currently do this. Failing to estimate partition functions and strength 
functions consistently can lead to inaccurate predictions of hnal equilibrium parameters. In 
equilibrium, the hard won rates no longer matter and only the partition functions govern 
the hnal quantities of interest {Ye, abundances, etc.). 


6. Results and Discussion 

6.1. A validity test: Comparison with Shell Model Based Rates for A=60-65 

Here we address the reliability of our calculated rates. We have shown in section 3 that 
with a simple prescription some gross features of the strength distribution, in particular the 
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total strength and centroid of the distribution, may be estimated. However, shell model and 
experimental results typically show a rich structure in the strength distribution, with this 
structure varying markedly from nucleus to nucleus (e.g. Caurier et ah (1999)). Do the rates 
depend sensitively on the hner features of the strength distribution? Or, alternatively, can 
a computationally simple method give a good estimate of the weak interaction rates? These 
questions can be addressed by comparing our relatively simply derived rates with rates based 
on more detailed large dimension shell model calculations. 

In hgures 7 and 8 we show the log of the ratio of our calculated beta decay rates to those 
calculated in Langanke and Martinez-Pinedo (2000) at several temperature/density points 
relevant for stellar collapse, and for all nuclei in the range A=60-65 for which LMP calculated 
rates. The horizontal axis in this hgure is the log of the beta decay rate calculated by LMP. 
We have presented the comparison in this way because nuclei with very small rates will 
not be so important in determining the evolution. The comparison is generally remarkably 
favorable, with typical results differing by less than a factor of two. Figures 9,10 and 11 are 
the same as the previous three hgures, but now electron capture rates are being compared. 
Again, the comparison is favorable, the rates as calculated in the two schemes being within 
a factor of three or so. This is surprising given the potential sensitivity of electron capture 
rates to the placement and the width of the Gamow-Teller distribution. 

The fact that a simple prescription and rough estimates of the strength do a reasonably 
good job in getting the rates relevant for stellar collapse is not to say that all rates are 
accurately determined at all temperatures and densities with a simple model. Some of 
our calculated rates deviate signihcantly from more detailed calculations. This typically 
occurs for /?“(/?+) decay when thermal population of the states in a collective Gamow-Teller 
resonance results in an exponential sensitivity, or for lepton capture when the maximum 
lepton energies are just on the edge of being able to reach the resonance in the daughter. 
Typically this exponential sensitivity is accompanied by a very small rate, so that it is not 
very important for stellar evolution. 


6.2. Examples of some rates important in X-ray burst environments 

X-ray bursts arise from the thermonuclear burning of hydrogen and helium accreted onto 
the surface of a neutron star in a binary system (see Wallace and Woosley (1981)). Gharac- 
teristic temperatures during burning are a few hundred keV, and characteristic densities are 
pYe ~ 10®g/cm^. The creation of nuclei heavier than A 40 occurs via the rp-process, in 
which nuclei undergo (p,7) reactions until they approach the proton drip line and/or /?+/ec 
decay intervenes. The time for burning along the rp-process path is set to some extent by 
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the j3^/ec lifetime of a few important waiting point nnclei. Here we present a few examples 
of important rates. 

Three important waiting point nnclei with A<80 are the even-even (proton bonnd) 
nnclei ®^Ge, ®®Se, and ^^Kr (Schatz et ah 2001). For each of these nnclei the single proton 
captnre danghter (Z-|-1,N) is nnbonnd, while the next heaviest isotone (Z-|-2,N) is bonnd. 
At low temperatnres (Tg < 1.5) the weak rate most important for determining flow towards 
the valley of (3 stability is A^+(Z,N). For higher temperatnres an eqnilibrinm between (Z,N) 
and (Z-|-2,N) is reached, so that \(^+{Z + 2,N) is also important. As a typical example we 
discuss 3 gKr and it’s two-proton capture daughter ggSr. 

The first 2~^ excited state of ^^Kr lies at ~ 700 keV and is not significantly populated 
for temperatures Tg < 2.5. The ground state lifetime of ^^Kr is experimentally-determined. 
The low lying strength distribution in the daughter (^^Br) has also been measured. With the 
ground state lifetime and low lying strength distribution measured, the only missing piece of 
information is the strength distribution at excitation energies too high to be experimentally 
observed. Our rough estimate (Eq. 21) places a sizable portion of the resonance strength 
within the Q-value window for the decay. In order not to conflict with the experimentally- 
determined lifetime we push this strength up to Q=0 in our calculation. 

To quantify the uncertainty in the rates arising from the high lying resonance strength 
we plot in Fig. 12 the total rate + ec) for ’^^Kr. Also included in this hgure is the fraction 
of the total rate coming from Q=0 and above. The hgure shows that the experimentally- 
determined lifetime is sufficient for a determination of the thermal decay rate at the 10% level 
for pYe < lO^g/cm^. At pY^ = lO'^g/cm^, our simple estimate shows the high lying strength 
accounting for ~ 20% of the total weak rate. Because our method puts essentially all of the 
strength at Q=0, and very little above Q=0, it is unlikely that we have underestimated the 
contribution to the rate from the high lying resonance strength. 

The decay of ^"^Sr is more easily calculated because the ground state (or a low lying 
excited state) of the odd-odd N=Z daughter C^^Rb) is the IAS of the ground state of '^"^Sr. 
The large matrix element and Q-value (~ 10 MeV) for the Fermi decay mean that electron 
capture cannot compete with j3^ decay in x-ray burst conditions for this case. It is difficult 
to reliably estimate the contribution to the j3^ rate from GT"*" transitions, but a reasonable 
estimate is that the diffuse GT+ strength only decreases the ground state lifetime by at most 
10 - 20 %. 

By mirror symmetry, a low lying thermally populated state of ^^Rb corresponds to the 
IAS of the next even-even j3^ nucleus (’^"^Kr), so that the chain ^"^Sr —Rb —Kr is 
fast and dominated by Fermi transitions. The decay of the odd-odd N=Z nucleus (^^Rb 
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in this example) should is calculated as the decay of a thermally populated back resonance 
as discussed above. This is because only those parent states with an IAS in the even-even 
daughter nucleus decay via the Fermi transition. Since the partition function of the odd-odd 
parent increases more rapidly with temperature than the partition function of the even- 
even daughter, the decay rate decreases rapidly with increasing temperature. Another set 
of nuclei important for the rp-process with decays dominated by the Fermi transitions are 
those nuclei with (Z = N -|- 1, N = even). These have an IAS near or at the ground state of 
the daughter. 

Generally, for a given ground state lifetime, the electron capture rate is small and 
insensitive to the placement of strength within the Q-value window as long as the initial 
state electron energy is small compared to the energy of the decay, i.e. 

(A'e) < (Mp — Md -|- Ej — Ef). (25) 

In the x-ray burst environment {Ee) < 0.5 MeV, while the Q-values for the proton-rich nuclei 
are typically greater than 4 MeV, so the thermal lifetime is reliably estimated from the ground 
state lifetime. By contrast, nuclei closer to the valley of stability have smaller Q-values and 
therefore the lifetimes are nearly entirely determined by thermal electron capture. 

For example, consider ®®Ge, a nucleus important in steady state burning on and/or near 
neutron star surfaces. This nucleus has a Q-value of only ~2 MeV and a fair portion of 
the strength lies at the upper end of the Q-value window. Gonsequently thermal electron 
capture dominates over positron decay for pi/. > lO^g/cm^. This is shown in figure 13. 


6.3. Weak rates in the late time pre-supernova star 

Here we discuss some of the systematics of the weak rates in the hot and electron- 
degenerate core during the ~ lO^s before core bounce. Most of the nuclei with A>65 present 
in the core will be blocked or nearly blocked to GT^ transitions in the zeroth order shell 
model picture. Our calculation of the electron capture rates for these nuclei is based on an 
estimate of the GT^ configuration mixing strength. Forbidden transitions and transitions 
allowed by the thermal unblocking of strength will also contribute to the electron capture 
rates. It is useful to estimate how large the configuration mixing strength must be in order 
to justify the neglect of these latter types of transitions. 

Thermal unblocking refers to the population of parent excited states that have zeroth 
order shell model configurations consistent with allowed transitions to states in the daughter. 
Assume that these parent states comprise a fraction 5Z of the total partition function. 
Then, an effective thermal unblocking matrix element is |MtuP ~ /Z)\Ms,p\'^. Here 
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iMi.p.P ~ 1 — 3 is a typical single particle allowed transition matrix element. An accurate 
estimate of 5Z is very difficult and an important open issue. Fuller 1982 parametrizes 
5Z K, Zexp{—E*/T), with E* the excitation energy of the lowest parent excited state with 
an allowed GT+ transition in the zeroth order shell model picture. With this schematic 
notation, thermal unblocking can compete with a configuration mixing strength of 0.1 — 1 
if E*/T < 2. For a typical E* ~ 5MeV, then, thermal unblocking can be neglected for 
T < 2 MeV. for T > 2 MeV the thermal population of the GT~ resonance states (which we 
do include in our calculations) also becomes important. 

Forbidden transitions become important as the electron chemical potential increases 
and the wavelength of the leptons involved in an electron capture event become small 
enough to probe structure in the nucleus. Again following the convention of Fuller 1982, 
the contribution to the electron capture rate from forbidden transitions can be written as 
Afor ~ |Affor| /for ( A'for ) Q) P-e)• Here |iVffor| 10 — 20 is roughly the number of protons in the 
fp shell multiplied by a typical single particle hrst forbidden matrix element. The unique 
hrst forbidden phase space factor /for depends on the centroid in energy of the forbidden 
strength distribution Efov, the parent-daughter mass difference Q, and the electron chemical 
potential pe- For a typical Q = 10 MeV, forbidden transitions compete with a low lying 
conhguration mixing strength of ~ 1/2 for pe = 31 MeV if F^for = 5 MeV. If the centroid 
of the forbidden strength lies instead at ~ 10 MeV above the daughter ground state, then 
forbidden transitions do not become important until pe > 37 MeV. 

With the assumption that the high density electron capture rates are dominated by 
transitions involving low-lying conhguration mixing strength, these rates are trivial functions 
of the electron Fermi energy and the parent-daughter mass difference. This is shown in hgure 
14 where the electron capture rates are presented for all nuclei with (N-Z)/A>0.1 that fall 
in the mass range A=66-80. The dependence of the rate on the Q-value for the transition is 
given by the simple analytic expression 


\^2 ln(2)10 ^•®|MGTp(F4(?7efr) + g^F2(77eff) + 2gF3(?7efr)). (26) 

Here pefr = {Up + me + q — i7res,d)/F, with E^es,d the centroid of the GT+ resonance in the 
daughter, q = \q — i7res,d|/F, and the extra prefactor of 2 approximately accounts for the 
Coulomb distortion of the incoming electron for nuclei with Z ^ 30 — 40. This equation 
assumes that q — i7res,d is negative. The functions 


Fk(rrt 


7 °° x^dx 

Jo + 1 


(27) 


appearing in Eq. 26 are the relativistic fermi integrals. Setting E^es.d = 2 MeV and 
lAfcTp = 1 in Eq. 26 gives a reasonable estimate of the electron capture rates (when 
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the rates are appreciable) when U-p > 10 MeV. These expressions can be used in place of 
our rate tables under these conditions. Using the simple approximations to the fermi func¬ 
tions developed in Fuller, Fowler and Newman (1985) gives an analytic approximation to the 
high density electron capture rates that is accurate to within the uncertainty in these rates. 
Electron capture rates for neutron rich nuclei at high electron fermi energies are essentially 
a function of only 1 parameter, the total GT^ strength within a few MeV of the daughter 
state. Discrepancies between our calculated rates and more detailed calculations, e.g. monte 
carlo-|-RPA calculations (Langanke, Kolbe, & Dean 2001), reflect the difference between our 
adopted |MgtP = 1 and the more detailed estimate of Mqt- 

Fig. 14 also gives an estimate of the uncertainties in the electron capture rates. At 
Tg = 10, pYf, = 3 ■ 10^°g/cm^ (Uf ~ 16MeV), an error of 2 MeV in the position of the 
centroid of the strength results in a change in the rate by a factor of ~ 4 for the nuclei with 
the smallest rates, and a factor of ~ 2 for the nuclei with the largest rates. At Tg = 10, 
pYf. = lO^^g/cm^ (Uf ~ 23 MeV), the uncertainty in the placement of the centroid implies 
an error of at most a factor of two in the electron capture rate. The uncertainty in the total 
strength is probably about a factor of 2 on average. 

The systematics of the (3 decay rates at high temperature and density are also simple. 
In Fig. 15 we show the (3 decay rates at Tg = 10, pYe = 10^°g/cm^ for all nuclei with 
(N — Z/A) > 0.1. These rates fall into three distinct bands corresponding to odd-odd nuclei, 
even-even nuclei, and even-odd/odd-even nuclei. This can be understood by noting that at 
Tg = 10 the GT+ resonance is thermally populated (under our assumption that the strength 
is centered at 2 MeV for these blocked nuclei). In this case the rate of decay is approximately 
A(Q)Zdaughter/-^parent, where A(Q) has a simple dependence on the Q-value. For even-even 
parents, 2’daughter/^parent U typically approximately 20 at Tg = 10, for odd-odd parents the 
ratio is about 1/20, and for odd-even/even-odd parents the ratio of partition functions is 
about 1. Note that the f3 decay rates are more sensitive to the (unknown) details of the 
strength distribution than the electron capture rates. At Tg = 10, an error in the placement 
of the centroid of the resonance of 2 MeV changes the beta decay rate by about an order of 
magnitude. 


6.4. Discussion and Conclusions 

We have provided estimates of weak interaction rates for nuclei in the mass range A=65- 
80. These may be useful in simulations of x-ray bursts and pre-supernova stellar evolution. 
The rates have been calculated using available experimental information and simple estimates 
for the strength distributions and matrix elements for allowed and discrete-state transitions. 
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The efficacy of our approach is conhrmed through comparisons with detailed shell model 
based rates for nuclei in the mass range 60-65. 

The single most uncertain aspect of the rate calculations is the GT"*" resonance strength 
distribution. Our simple prescription gives good overall agreement with detailed shell model 
calculations of the GT"*" strength for nuclei with A<65. However, our prescription is probably 
not reliable for nuclei at the end of the fp shell. Fortunately, nature to some extent does 
not seem to care about some of the hard to get details of the strength distribution for these 
nuclei. In the pre-collapse supernova this is because the nuclei present are nearly blocked to 
GT^ transitions. (n,p) exchange experiments on such nuclei show that the strength is broad 
and low lying in the daughter. Because the electron fermi energies are high in the dense 
pre-supernova Fe-core, this implies that the electron capture rate is principally a function 
of the experimentally-determined parent-daughter mass difference. However, our work does 
not provide the detailed estimates of the magnitude of the configuration mixing strength 
that will ultimately also be needed for supernova simulations. The (3 decay rates for the 
nuclei present in the late-time pre-supernova core are exponentially sensitive to the centroid 
energy of the GT"*" resonance and are less certain. In x-ray burst environments the opposite 
condition holds: the parent daughter mass differences are typically large compared to the 
electron energies, so that the decay rate is dominated by the experimentally-determined 
lifetime. 

The proper treatment of the high temperature electron capture and particulary beta 
decay rates is a challenging and important issue. Special care is needed in evaluating the 
partition functions for these rates. The beta decay rates for nuclei in these conditions is 
determined, among other things, by the nuclear partition function. In turn, the equilibrium 
nuclear composition is determined by the competition between beta decay and electron cap¬ 
ture. If the partition functions used to estimate the composition at a given Y^. do not match 
with the partition functions used to calculate the weak rates, the calculated equilibrium of 
the sytem will not be correct. 
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Table 1. Comparison of the shell model GT“ centroids calculated in Langanke (LMP) 
with those estimated from Eq. 20 (Present). All energies are in MeV and refer to energy 

with respect to the daughter ground state. 


Parent Nucleus 

LMP 

Present 

55 Fe 

12.6 

10.4 

56pga 

9.6 

8.5 

STpe 

12.6 

11.5 

58Fe 

11.0 

9.5 

59 Fe 

13.6 

12.6 

60 Fe 

10.3 

10.6 

61 Fe 

13.8 

13.75 

62Fe 

11.8 

11.7 

58f^ja 

9.2 

6.44 

59Ni 

10.6 

9.46 


9.0 

7.5 

61 Ni 

13.3 

10.6 

62 Ni 

9.2 

8.2 

63 Ni 

13.2 

11.5 

64 Ni 

9.6 

9.15 

65 Ni 

12. 

12.35 

66 Co 

13.2 

12.7 

6iCo 

12.5 

10.77 

68 Co 

14.7 

13.7 

69 Co 

13.1 

11.6 

60 Co 

14.0 

14.8 

61 Co 

13.6 

12.6 

62 Co 

15.3 

13.9 

63 Co 

14.4 

13.7 

64 Co 

16.0 

16.7 

66 Co 

14.6 

14.8 

66Mn 

13.0 

11.8 

66Mn 

14.7 

15.1 

61'Mn 

13.1 

13.5 

68Mn 

15.5 

16.4 

59Mn 

14.0 

14.7 

60 Mn 

16.0 

17.4 

6lMn 

16.2 

15.5 


^For these nuclei the results from 
(p,n) experiments give a centroid ap¬ 
proximately 0.7-1 MeV lower than the 
shell model results. 


Note. — The shell model centroids in 
this table were taken from the tables or 
estimated from the graphs in Caurier et 
al. 1999 and LMP. 
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Table 2. Comparison of the shell model GT"*" centroids calculated by LMP with those 
estimated from Eq. 21 (Present). All energies are in MeV and refer to energy with respect 

to the daughter ground state. 


Parent Nucleus 

LMP 

Present 

55 Mn 

4.6 

4.6 

56Mn 

5.9 

6.4 

56 Fe 

2.6 

2.4 

56 Co 

8.2 

00 

bo 

58Mn 

5.5 

5.15 

58 Co 

7.35 

8.1 

58Ni 

3.75 

3.65 

59 Co 

5.05 

5.0 

60 Co 

6.35 

6.74 

60 Ni 

3.4 

2.7 

61 Fe 

2.1 

1.8 

61 Co 

3.7 

3.4 

6lNi 

4.7 

4.7 

6iCu 

6.7 

6.4 

62 Ni 

2.1 

1.8 

64 Ni 

1.3 

1.8 


Note. — The shell model centroids 
in this table were also taken from Cau- 
rier et al. 1999 and LMP. Where an ex¬ 
perimental result is listed along with the 
shell model result we have presented the 
experimental result. 
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Thermally populated states in the GT- resonance undergo 
fast transitions to low lying states in the daughter. 



Fig. 4.— Illustration of the position and energetics of discrete state-resonance transitions 
(back resonances). 
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Fig. 5.— Comparison of our decay rates with the LMP rates for nuclei in the mass 
range A=60-65. The comparison is made at Tg = 30, pYe = lOg/cm^. These conditions are 
artihcial but serve to illustrate the influence of the treatment of the partition function on 
estimates of the weak rates. In all plots where a rate is shown the rate is in units of sec“^. 
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Fig. 6.— Same as Fig. 5 but at Tg = 10, pYe = 10®g/cm^, conditions approximating the 
post Si burning degenerate core before collapse. 
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Fig. 7.— Comparison of our beta decay rates calculated using a simple estimate of the 
strength distribution with the more sophisticated calculations of Caurier et ah (1999) and 
Langanke and Martinez-Pinedo (2000). The different circles correspond to different nuclei. 
Comparisons for all nuclei in the mass range A=60 — 65 for which Langanke and Martinez- 
Pinedo (2000) provide rates are presented. Here Tg = 3 and pYe = lO^g/cm^. 
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Fig. 8.— Same as figure 7 but with Tg = 5 and pY^ = 10®g/cm^ 
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Log(LMP rate) 

Fig. 9.— Comparison of our electron capture rates with those of Langanke and Martinez- 
Pinedo at Tg = 3, pY^ = lO^g/cm^. Only results for nuclei that are estimated by LMP to 
have electron capture rates larger than 10“^°sec“^ are shown. 
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^ig. 11.— Same as figure 9 at Tg = 10 and pY^ = 10®g/cm^. 
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Fig. 12.— Influence of high lying resonance strength on the total rate {j3^ + ec) for ^^Kr. In 
this hgure /high is an estimate of the fraction of the total rate arising from transitions involving 
experimentally unmeasured, high lying resonance strength. The upper curve, labelled A/Aq, 
shows the ratio of the total rate to the j3^ decay rate at zero temperature. 
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Fig. 13.— /ec rates for ®®Ge.The increase of the j3^ rate with temperature arises from the 
fast decay of the thermally populated hrst .P = 2+ excited state. Unlike the O’*" ground state, 
the 2+ state can have allowed transitions to several 2+ and d"*" daughter states. The decrease 
in the electron capture rate at low temperatures is likely an artifact of our calculation and 
arises because the O’*" ground state has substantial GT strength at high daughter excitation 
energies. As the daughter level structure is poorly known at these high excitation energies, 
our calculation does not include estimates of the allowed matrix elements for the decay of 
the hrst excited 2+ state to those states. 
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Fig. 14.— Electron capture rates for nuclei in the mass range A< 65 < 80 and for which 
(N-Z)/A>0.1 as a function of the Q-value (dehned here as the parent daughter atomic mass 
difference). These rates have a simple dependence on the Q-value (see Eq. 26) because of our 
assumption that the conhguration mixing strength is nucleus-independent for GT^ blocked 
nuclei. 
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Fig. 15.— f3 decay rates for neutron rich nuclei in the mass range A< 65 < 80 at Tg = 10 
and pYe = 10^°g/cm^. The systematics illustrated here are discussed in the text. 




